The family of strange multiquarks related to the D s (2317) and D s (2457) 
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I study the IX, (2317) and D s (2457) discovered at BABAR, CLEO and BELLE, and find that they 
belong to a class of strange S = — 1 tetraquarks and pentaquarks, which is equivalent to the class of 
kaonic molecules bound by short range attraction. In this class of hadrons a kaon is strongly trapped 
by a s-wave meson or baryon. To describe this class of multiquarks the Resonating Group Method 
is applied to a standard quark model with chiral symmetry breaking, and the short range kaon- 
meson(baryon) interactions are extracted. A criterion is derived to classify the attractive channels. 

I conclude that the mesons Bi° + \ B^ 1 ^ , and the baryons S7 CC , fici, &bb clearly belong to the new 
hadronic class of the Z) s (2317) and D 3 (2457). The hadrons /o(980), A, E c , E;, possibly belong to a 
related family. 



INTRODUCTION 



Recently new narrow scalar resonances D s (23171 and 
D s (2457) were discovered at BABAR Q, Cleo and 
Belle 0. They are not expected to be standard hadrons 
because their masses are close to hundred MeV smaller 
than the ones predicted in quark models for positive par- 
ity D s quark-antiquark mesons For instance the ex- 
perimental candidate to the first quark-antiquark D s me- 
son is the An (2536) @. This is confirmed by lattice 
simulations 5j for quark-antiquark mesons. The exis- 
tence of non-mesonic and non-hadronic multiquarks has 
been suggested form the onset of the quark model 0, H| • 
These positive parity £>,, resonances are interpreted by 
Beveren and Rupp [jj as poles in the S Matrix of the 
coupled channels of mesons and meson pairs. They were 
also predicted by Nowak, Rho and Zahed [TollTlT ] as chiral 
partners of the well known negative parity pseudoscalar 
D s (1968) and Z? s (2112). They are as well interpreted as 
K-D molecules by Barnes, Close and Lipki n 1131 . simi- 
lar to the Deusons anticipated by Tornqvist In the 
same way I find here that the new D s resonances can 
be understood as tetraquarks, or equivalently as s-wave 
D-K molecules bound by the short range attraction. 

The experimental discovery of these hadron also sug- 
gests the existence of a new class of multiquarks. In this 
paper I study the family of all possible narrow tetraquark 
and pentaquark resonances where the quark s, or the 
S = — 1 Kaon play a crucial role. The strangeness 5=1 
pentaquark + is more difficult to bind and was recently 
addressed with the same techniques of this paper in ref- 
erence 14] . A chiral invariant framework fla , ITtI ITsj 
is used to compute microscopically, at the quark level, 
the masses of this new class of hadrons. The resulting 
mechanism which provides the binding of the S = — 1 
tetraquarks and pentaquarks is equivalent to the short 
range attraction of hadrons. 

Here a standard Quark Model (QM) Hamiltonian is 



assumed, 



(i) 



i<j 



where each quark or antiquark has a kinetic energy Ti 
with a constituent quark mass, and the colour depen- 
dent two-body interaction Vij includes the standard QM 
confining and hypcrfine terms, 
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-A; • A,- 



V con f(r) + V hyp (r)S l ■ Sj 



(2) 



Moreover the Hamiltonian includes a quark-antiquark 
annihilation term which is the result of spontaneous 
chiral symmetry breaking. 

This paper is organised in sections. In Section II the 
QM is reviewed, together with the Resonating Group 
Method (RGM) [19( which is adequate to study multi- 
quark states, where several quarks overlap. The RGM, 
together with chiral symmetry, produces short range 
hadron-hadron potentials, which can be either repul- 
sive (hardr core repulsion) or attractive. In Section 
III a criterion is derived to discriminate which sys- 
tems bind and which are unbound. This criterion is 
applied to find, among the s-wave hadrons, the can- 
didates to trap a kaon. In section IV the binding 
energy is computed for the selected positive parity 
mesons /o (980) (°+),^ 0+) , D { s 1+ \ £^ 0+) , B { s 1+) and neg- 
ative parity baryons A, S c , S5, Sl cc , f2 c &, £lbb- Finally the 
results are presented and discussed in Section V. 



II. STUDYING MULTIQUARKS WITH THE 
RGM 



For the purpose of this paper the details of potential 
are unimportant, only its matrix elements matter. 
The hadron spectrum constrains the hyperfine potential, 



hyp) 



l(M A 



M N ) 



(3) 
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The quark-antiquark annihilation potential Aq is also 
constrained when the quark model produces spontaneous 



FIG. 1: The Jacobi coordinates of the incoming four quark 
wave-function 4>a(pa) <\>b{pb) iP(\a,b)- 



chiral symmetry breaking |l5l Il6j | . The annihilation po- 
tential A is present in the 7r Salpeter equation, 



2T + V A 
A 2T + V 



(4) 



where the tt is the only hadron with a large negative en- 
ergy wave-function, §~ ~ 4> + . In eq. Q the annihilation 
potential A cancels most of the kinetic energy and confin- 
ing potential 2T + V. This is the reason why the pion has 
a very small mass. From the hadron spectrum and using 
eq. Q the matrix elements of the annihilation potential 
are determined, 



(2T + V) s =v 
=> {A)s=o 



-(2M N -M A ) 
-~(2M N -M A ) 



(5) 



where this result is correct for the annihilation of u or 
d quarks. When a strange quark is present, the corre- 
sponding matrix element is smaller by a factor a which 
is a power of the constituent quark mass ratio M Ui d/M s 
The QM of eq. (IJ reproduces the meson and baryon 
spectrum with quark and antiquark bound-states (from 
the heavy quarkonium to the light pion mass). Moreover 
the RGM was first used in hadronic physics by Ribeiro 
pof to show that in exotic hadron-hadron scattering, the 
quark-quark potential together with the Pauli repulsion 
of quarks produces a repulsive short range interaction. 
For instance this explains the N — N hard core repul- 
sion, preventing the collapse of nuclei. This N — N 
hard core repulsion is supposed to also occur in several 
hadron-hadron interactions, and this explains why many 
multiquarks systems are not stable. However Deus and 
Ribeiro [2l| used the same RGM to show that, in non- 
exotic channels, the quark-antiquark annihilation could 
produce a short core attraction. Recently it was shown 
that in the particular case of the low energy tt — t: system 
in the chiral limit, the short range attraction and repul- 
sion exactly cancel |17| , resulting in a Adler Zero and the 
Weinberg theorem. Addressing a tetraquark system with 
the 7T — 7r quantum numbers, it was shown that the QM 
also fully complies with the chiral symmetry, including 
the PC AC theorems ^?|]- Therefore the QM is adequate 
to address the new D s hadrons, which were predicted by 
Nowak, Rho and Zahed in an effective chiral model. In 
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FIG. 2: Examples of RGM overlaps are depicted, in (a) the 
norm overlap for the meson-baryon interaction, in (b) a ki- 
netic overlap the meson-meson interaction, in (c) an interac- 
tion overlap the meson-meson interaction, in (d) the annihi- 
lation overlap for the meson-baryon interaction. 



this paper the QM and the RGM are applied to S = — 1 
multiquarks. 

The RGM computes the effective multiquark 

energy using the matrix elements of the microscopic 
quark-quark interactions. Any multiquark state can 
be decomposed in combinations of simpler colour sin- 
glets, the baryons and mesons. This can be illustrated 
with the colour structure of a multiquark with two 
quark-antiquark pairs, qi, q2, q%, q^. Assuming that this 
tetraquark is a colour singlet, each quark-antiquark pair 
can either be a colour singlet or a colour octet. Nev- 
ertheless the octet-octet state can be described with an 
exchange operator and with the singlet-singlet state, 



'1,2 ' °3,4 



27=2 ( Pl3 " V 



ll,2 1 3 ,4 , (6) 



and there is only one anti-symmetrised state, 

(1 -Pl 3 )(l - P 24 ) ll,2 1 3 ,4 = 



-J=(l-P 13 )(l-P24)8l,2 



53,4 



(7) 



where (1 — Pi3)(l — P24) is the quark and antiquark anti- 
symmetriser. Therefore the multiquarks can be described 
with an anti-symmetrised basis of baryons and mesons. 
The multiquark may only be a bound state, or a narrow 
resonance, if these baryons and mesons are sufficiently 
attracted. Thus the problem of multiquark stability can 
be technically reduced to the problem of the binding of 
baryons or mesons. 

The RGM produces both the energy of a multiquark 
state and the effective hadron-hadron interaction. The 
energy of the multiquark is computed with the matrix 
elements of the hamiltonian (JIJ. The wave functions of 
quarks are arranged in anti-symmetrised overlaps of sim- 
ple colour singlet hadrons. As an example, in Fig. ^ 
a tetraquark system is arranged in a pair of mesons A 
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and B. Once the internal energies E A and Eb of the two 
hadronic clusters are accounted, 



= E a +E b + V at 



(8) 



where ^2 p (~^) p P is the anti-symmetrizer, the remain- 
ing energy of the meson-baryon or meson-meson system 
is computed with the overlap of the inter-cluster micro- 
scopic potentials, 

V* Z A = {4> B (p A \-{VlA + Vi 5 + 2V 2 4 + 2V 25 )3Pl4 

+3A 16 \4> a <I>b)/{<I>b M 1 ~ ZPu\4>a4>b) 
V^sA = {<Pb(I>a\(1 + Pab)[-(Vi3 + V23 + Vu + V2±) 

mcs B 

xPl3 + A 23 + A U ]\(j) A <pB) 

/{4> B (j> A \{l + Pab){1 - Pi 3 )\Mb) , (9) 

where Pij stands for the exchange of particle i with par- 
ticle j, see Fig. 



III. A CRITERION FOR BINDING 

In what concerns the annihilation potential, it only 
occurs in non-exotic channels. Then it is clear from eqs. 
(|5|) that the annihilation potential provides an attractive 
(negative) overlap. This confirms that the hard core can 
be attractive for non-exotic channels where annihilation 
occurs. 

In what concerns the quark-quark(antiquark) poten- 
tial, it may also contribute to exotic channels. Because 
the potential Vij is assumed to be proportional to the 

colour dependent 4f • -4- it is clear that it can only con- 
tribute together with an exchange interaction, which pro- 
vides a color octet, see eq. Moreover the spin inde- 
pendent part of the interaction vanishes. For instance 
in the meson-meson overlap of eq. © the overlap of 
2 ' 2 2 ' 2 essentially cancels with the overlap of 
2 "2 + 2 "2 • only potential which may contribute 

is the hyperfine potential, proportional to %■ • S{ ■ Sj . 
In the present case where the kaon is a spin singlet, the 
minus phase from colour is inverted by a minus phase 
from spin. I find that the total colour and spin matrix 
element is a hyperfine splitting, 

(^13(^13 + V 23 + V u + V 24 )) = 4^(M A - Mjv) . (10) 

The flavour trace is quite simple and the spatial integral 
provides a geometrical overlap. All the corresponding 
factors are positive, an therefore the quark-quark inter- 
action results in a repulsive interaction. 

These results are independent of the particular quark 
model that one chooses to consider, providing it is chi- 
ral invariant. Therefore two opposite classes of diagrams 



exist. The exchange diagrams produce a repulsive inter- 
action, which turns out to be proportional to the hyper- 
fine quark-quark(antiquark) interaction. The annihila- 
tion diagrams produce an attractive interaction, which 
turns out to be proportional to the spin-independent 
quark-quark(antiquark) interaction. I arrive at the at- 
traction/repulsion criterion for the short range hadron- 
hadron interaction, 

- whenever the two interacting hadrons have quarks (or 
antiquarks) with a common flavour, the repulsion is in- 
creased, 

- when the two interacting hadrons have a quark and an 
antiquark with the same flavour, the attraction is en- 
hanced . 

This paper is dedicated to the class of resonances which 
can be understood as a S—-1 kaon su or sd strongly 
trapped by a s-wave hadron. The criterion shows that 
all hadrons with an antiquark u or d or with a a quark 
s allow the exchange overlap (Pi 3) , and this would cer- 
tainly contribute to repulsion [23.1231 124| . Although sys- 
tems with both a short range repulsion and a short range 
attraction may still bind |14| , I specialize here in systems 
where clearly there is no hard core repulsion. In what 
concerns attraction, a quark u or d is needed in the s- 
wave partner of the kaon, in order to produce annihi- 
lation. Therefore the isospin of the partner of the kaon 
needs to be close to the opposite isospin of the kaon. This 
excludes the mesons rj, rf ui , (/)■■■. Moreover I also spe- 
cialise in systems which may achieve a remarkable sta- 
bility, where the kaon partner is a hadronic resonance 
with a very narrow width. This not only essentially re- 
stricts the kaon partner to s-wave mesons and baryons, 
it also excludes the meson p and the baryon A because 
they are very wide, due to the decay with a pion produc- 
tion. The pion is also excluded because it is too light to 
bind to the kaon, all that one may get is a very broad 
resonance, the kappa resonance |25|. which has been re- 
cently confirmed by the scientific community. There- 
fore the hadrons which are best candidates to strongly 
bind the Kaons si are the s-wave hadrons with flavor 
Is, Ic, lb, III, lie, lib, lec, leb, Ibb. This is expected to re- 
sult in the / (980) the £> s (2320 @ the D s (2463) Q, 
the A(1405) and several other resonances. 



IV. THE BINDING ENERGY 

A convenient approximation for the meson and baryon 
s-wave wave-functions is the harmonic oscillator gound- 
state wave-function, 



^000 ip P ) =Nct : exp 



Vp 
2a 2 



a 



(11) 

where, in the case of vanishing external momenta pa and 
Pb, the momentum integral in eq. Q is simply Af a ~ 2 - 
The coordinates of the incoming 4>a<Pb functions are il- 
lustrated in Fig. ^ while the coordinates of the outgoing 
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E 



turn out to be multiplied by the gaussian separable fac- 




FIG. 3: The matrix element of the Green function go(E, [i, a) 
are plotted in dimentionless units of /i = a = 1. The general 
case scales as go(E.y) = -^gl{-^E, 1). 



4>\<l>\t have the quark 1 and 3 exchanged. I summarise 
O EE H3 the effective potentials computed for the dif- 
ferent channels, 



Vis--* 


= (2 + a)\{A)Af a 

6 




= 2^{A)N a - 2 , 




= 2^{A)N a - 2 , 




= 2^{A)N a - 2 , 


Vr-b* 


= 2 l -{A)N a - 2 , 
6 


Vk-n 


= Z-(A)N a - 2 , 
6 


Vk-h € 








Vk-3 cc 


6 


V K -B cb 


6 


V K -s bb 


= 2 l -(A)N a - 2 , 
6 



(12) 



where the colour and spin factors contribute respectively 
with 1/3 and 1/2, (A) is of the order of 430 MeV and 
the geometrical factor is Af a ~ • The remaining factor 
is the flavour factor. The parameter a is the only one 
that is model dependent, and it will be determined with 
a fit of experimental binding energies. The estimation 
of a is an important by-product of this method because 
the hadronic size can not be estimated directly by the 
hadronic charge radius which is masked by the vector 
meson dominance. 

To study binding one has to proceed to the finite mo- 
mentum case. Then the effective potentials in eq. I|12f) 



tor, exp 



Pa 
~2/3 2 



f rf3 PA 

J (2tt) 3 



exp 



Pa 
" 2/3 2 



In the exotic chan- 



nels with exchange diagrams only, this result can be 
prooved j2||, and moreover the new parameter j3 = a. 
This occurs because the overlaps decrease when the rel- 
ative momentum of the hadrons A and B increases. In 
the non-exotic channels with annihilation diagrams the 
present state of the art of the RGM does not allow a 
precise determination of the finite momentum overlap. 
Nevertheless, it is expected that eventually the overlap 
decreases due to the geometrical wave-function overlap 
in momentum space. Here the precise determination of 
f3 does not affect the results. I order to reduce the num- 
ber of parameters, and for simplicity, in this paper I as- 
sume that j3 ~ a. For other approaches that also lead 
to a separable potential, see for instance Ref. 0. This 
parametrisation of the Schrodinger equation in a separa- 
ble potential, 

[(E-T A -T B ){l + n\r ><4> a \) 
+ W |0 Q X <H] >= , (13) 

enables the use of standard techniques ^3 to exactly 
compute the scattering T matrix, 



T = \(/> a > 



1 



-nSo 



< 



(14) 



The binding occurs when the T matrix has a pole for a 
negative relative energy, 



1 - 



1-71 



gS(E,n)=0. 



(15) 



In Fig. [3] the function Go is plotted. Go is real for a 
negative relative energy E and is complex for positive E. 
The binding only occurs if, 



(16) 



Using Fig. [3] it is possible to determine the parameter a 
which reproduces the experimentally measured binding 
energy of the D s (2320). With a binding energy of 46 
MeV, the corresponding parameter a is 285 MeV. This 
corresponds to a radius of 0.7 Fm for the nucleon. 



V. RESULTS AND CONCLUSION 

I now compute the binding energies of the kaonic 
strongly bound molecules, or equivalently strange multi- 
quarks. These multiquarks are divided into two different 
families which are respectively coupled and decoupled to 
pionic channels. The computations are straightforward 
and the results are displayed in Tables HI and ITT1 
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channel 



V2 



D s (2317) = 
D s (2457) = ■ 

g _ K~B" ! h " D 

B s = 

n+ = 



>/2 



.V2 



n2» = 



%/2 



Mea;p ^th Q! — /3 -Btfr -Beajp 



392 


-143 


285 


4 u 


46 


398 


-143 


285 
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46 


453 


-143 


285 


55 




454 


-143 


285 


55 




442 


-143 


285 


53 




466 


-143 


285 


56 




475 


-143 


285 


58 





channel 

_ K~K + +K"K" 



TABLE I: This table summarises the parameters /i, v ,a ,/3 
and binding energies B (in MeV) [6j for the channels closed 
to pion decay. The italic binding energy Bth of the D s (1327) 
is fitted from experiment. 



A. A family of very narrow multiquark resonances 



The dI° +>> has a lower mass by a hundred or more MeV 
than the expected 3 Po excitation of the Ds ''(1968). 
In order to conserve the angular momentum and par- 
ity, it can only decay to a Ds '(2112) with the cre- 
ation of a pion, in a p-wave, with a very low energy 
which is suppressed by an Adler zero 1171 1 . More im- 
portantly, this decay mode violates isospin conservation, 
and therefore it is quite suppressed. The same isola- 
tion from other strong hadronic channels is common to 
other multiquarks, with a quark s, one or more heavy 
quarks and the isospin zero combination uu + dd. I find 
that the Di 0+ \ L>i 1+) , Bf + \ B { s 1+) belong to the same 
class of tetraquark hadronic resonances. This class is 
equivalent to the picture of a kaon trapped by a s-wave 
meson with a short range attraction. In what concerns 
pentaquarks, where the heavy antiquark in the trapping 
meson is replaced by a pair of heavy quarks, the re- 
sults of this paper predict that there is a similar bind- 
ing with the quantum numbers of the fl cc , &cb, and 
flbb- I conclude that the tetraquarks and pentaquarks 
Di 0+ \ D { s 1+ \ Bi 0+ \ B { s 1+) ,n cc , n cb and fl bb belong to 
the same family of very narrow, non-exotic multiquarks. 
The corresponding binding energies of the kaon-hadron 
system are shown in Table [I] 



B. A second family, coupled to pionic channels 

A second family of tetraquarks and pentaquarks may 
exist, where the coupling to pionic channels does not vi- 
olate the conservation of isospin. The only suppression 
that one may expect comes from the Adler Zero. The 
only tetraquark that we select in this class is the 1 = 
fo, which is coupled to the tt — ir channel. A possible 
existing candidate is the /o(980), although it is not ex- 
cluded that it is a simple qq meson 0] . In what concerns 
pentaquarks, when the baryon (that traps the kaon) has 
two or more light quarks, the corresponding pentaquark 
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285 
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68 




456 


-167 
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75 




456 


-167 


285 
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TABLE II: This table summarises the parameters /i, v , a , j3 
and binding energies B (in MeV) 0] for the channels open to 
pion decay. 



is again coupled to pionic channels. For instance the kaon 
can be trapped by a nucleon to produce a pentaquark A, 
and this is coupled to the s-wave S s — tt channel, with 
the same isospin. A possible candidate to this state is 
the A(1405). Similarly the kaon can be trapped by the 
E c or the E& baryons to produce pentaquarks with the 
quantum numbers of negative parity H+, HJj, S°, Sj~~. 
These channels are respectively coupled to the tt — S c 
and 7r — Sf, channels where the S c and Hj have a positive 
parity. Nevertheless I compute the binding energies of 
this class of multiquarks, ignoring the effect of the cou- 
pled pionic channel. The results are displayed in TablelTTl 
I find that the binding energy may be excessively large 
in the fo and A channel. This suggests that there is mix- 
ing between the two coupled multiquark states, and that 
this pushes the mass of the fo and A up to the higher 
experimental results. For instance it been advocated by 
Oset and collaborators [28| that the A(1405) results from 
the mixing of a singlet-singlet state with an octet-octet 
state. 



Comparing multiquarks with molecules and 
with chiral excitations 



The formalism used here is also convenient to ad- 
dress the question, are these new hadrons multiquarks, 
or molecules composed of standard hadrons, or chiral im- 
ages of standard hadrons? 

I find that, in a microscopic and chiral invariant calcu- 
lation, the new hadrons appear as tetraquarks, or as pen- 
taquarks. In a macroscopic hadron-hadron interaction 
perspective, the resulting mechanism which provides the 
binding is the short range strong attraction of hadrons. 
Technically is is not possible distinguish a microscopic 
multiquark from a macroscopic strongly bound s-wave 
molecule, because the hadrons totaly overlap. The purely 
molecular pattern only appears in different cases, when 
there is some repulsion which produces the clustering of 
quarks, and this may happen in the "pentaquark" fam- 
ily of exotics 0|- Nevertheless the molecular perspec- 
tive is convenient to estimate the mass of the multiquark 
because the binding energies are not very large, and the 
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mass is essentially the sum of the kaon plus hadron mass, 
with a relatively small binding energy. 

The chiral excitation perspective, present in the Chi- 
ral Soliton Model ^(J or m the Chiral Lagrangian, p9| 
may also be eventually equivalent to the multiquark, or 
to the two-hadron narrow bound molecule. Suppose that 
a flavour singlet quark-antiquark pair uu + dd or ss is 
created in a given hadron H . When the resulting multi- 
quark H 1 remains bound, we have a state with an oppo- 
site parity to the original H. The reversed parity occurs 
due to the intrinsic parity of fermions and anti-fermions. 
Moreover, because s- waves are expected to have the low- 
est mass, this is equivalent to adding a pseudoscalar me- 
son to the original H . The finite masses of the Kaon 
and of the pion prevent the existence of a large num- 
ber of pseudoscalar mesons in the multiquark family of 
the -D s (2317). In this family the mass shift is of the or- 
der of the kaon mass minus the binding energy minus the 
strange- light quark mass difference — B — {M s — M u ), 
and this results in a mass shift close to 350 MeV. The 
same mass shift is expected in the whole family of the 
D s (2317). Thus the multiquarks H' studied in this pa- 
per can be regarded as the chiral partners of the original 
H, in agreement with Nowak, Rho and Zahed. Neverthe- 
less the standard 3 Po quark-antiquark D s meson remains 
a chiral partner of the 1 Sq one, where the mass shift is a 
result of the angular momentum and spin excitations. In 
the present framework both the quark-antiquark creation 
and the angular momentum excitation chiral excitations 
can be understood. 



D. Outlook 

The presented results only depend on the hadronic size 
a and not on the details of the quark-quark interactions, 



because the assumptions are quite simple. I also assumed 
that the inverse radius parameters a and (3 are identical 
for all channels, although small channel dependences are 
expected. Moreover I neglected the meson exchange in- 
teractions because they are expected to be smaller than 
the hard core interaction. In the same way I did not con- 
sider the s-channel coupling to a single meson or baryon. 
More importantly, in the second family of multiquarks, 
the coupling to pionic channels was neglected although 
this coupling conserves isospin. It would be particularly 
interesting to study these neglected effects in the second 
family studied in this paper, in order to check if they 
correct the computed binding which seems too strong. 
Because the research in this direction will be model de- 
pendent, this will be done elsewhere. 



The short range attraction studied in this paper is 
quite general and can also be applied to the different 
candidates to multiquarks which have been discovered 
quite recently |M El IH H HI E| . This will also 
be applied elsewhere. 
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